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(— I ' Abstract. We investigate the growth and structure of magnetic fields amplified by 

<*S |. kinematic dynamo action in turbulence with non-zero kinetic helicity. We assume a 

q | simple Gaussian velocity correlation tensor, which allows us to consider very large 

i—{ • magnetic Reynolds numbers, up to 10 12 . We use the kinematic Kazantsev-Kraichnan 

model of dynamo and find a complete numerical solution for the correlation functions 
of growing magnetic fields. 
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The prevailing theory for the origin of cosmic magnetic fields in the Universe is magnetic 
dynamo action, which is random stretching of magnetic field lines by motion of highly 
conducting plasmas or fluids in which these lines are frozen (e.g., Brandenburg and 
Subramanian 2005; Kulsrud 2005; Kulsrud and Zweibel 2008; Lynden-Bell 1994; Parker 
1979; Vainshtein and Zeldovich 1972; Zweibel and Heiles 1997). Magnetic fields in 
astrophysical systems are often observed to be correlated at scales significantly larger 
than the correlation scales of plasma motions. Such large-scale magnetic fields can be 
generated by dynamo action if the velocity field v(x, t) possesses nonzero kinetic helicity 
i:=/v-(Vxv)d 3 x^0 (Moffatt 1978; Steenbeck et al. 1966). 

In this work we explore the growth and structure of magnetic fields generated by 
dynamo action driven by a velocity field with non-zero kinetic helicity, in the framework 
of the kinematic Kazantsev-Kraichnan model. For simplicity, we assume a Gaussian 
correlation tensor of velocity fluctuations, which means that all velocity fluctuations are 
on a single-scale, that is, the turbulence consists of similar-size eddies. This assumption, 
although not commonly satisfied in astrophysical systems, allows us to consider very 
large magnetic Reynolds numbers, up to 10 12 . More general settings can be addressed 
by the same method, and will be investigated elsewhere. Fast dynamo action driven by 
Gaussian-correlated velocity field with non-zero helicity was first studied numerically in 
the framework of the Kazantsev-Kraichnan model by Kim and Hughes (1997). Their 
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work was limited to a study of magnetic field generation at very small scales, much 
smaller than the velocity correlation scale. In particular, they focused on the dependence 
of the fastest growth rate of the field on the value of the magnetic Reynolds number. 
Our work is complementary to their work. We consider magnetic field generation and 
find field growth rates at all scales. Thus, we present a complete numerical solution 
for helical dynamo action. In addition, a different numerical method that we use to 
solve the helical dynamo equations, allows us to consider magnetic Reynolds numbers 
much larger than whose investigated by Kim and Hughes (1997)fl In the next section 
we describe the Kazantsev-Kraichnan model. In the final section, we present our results 
and give our conclusions. 



2. Kazantsev-Kraichnan model of magnetic dynamo action 

We use the Kazantsev-Kraichnan model of kinematic dynamo action in homogeneous 
and isotropic turbulence (Kazantsev 1968; Kraichnan 1968). In this model, the ensemble 
statistics of velocity fluctuations is assumed to be Gaussian, with zero mean, (v) = 0, 
and the covariance tensor 

(?/(x, (x', £')> =« ij "(x - x')5(t - t'), (1) 

where is an isotropic tensor of turbulent diffusivity, 

613 -^)+ K ^+9^ k x k . (2) 

Here () denotes ensemble average, e* jfc is the unit anti-symmetric pseudo-tensor and 
summation over repeated indices is assumed. The first two terms at the right-hand side 
of ([2]) represent the mirror-symmetric, non- helical part, while function g(x) describes the 
helical part of the velocity fluctuations. For an incompressible velocity field (the only 
case we consider here), we have k^(x) = Kl(x) + xk' l (x)/2, where the prime denotes 
derivative with respect to x — |x|. Therefore, to describe the velocity field, we specify 
only two independent functions, kl(x) and g(x). 

The magnetic field correlator can similarly be expressed as 

S ij - — + Mi — — + Ke^ k x\ (3) 

X J X 

where the field solenoidality constraint divB = implies Mr?(x,t) = Ml(x, t) + 
(x/2)M' L (x,t). To fully describe the magnetic correlator, we therefore need to find only 
two functions, Mi(x,t) and K(x,t), corresponding to magnetic energy and magnetic 
helicity. The Fourier transformed version of equation ([3]) is 

<B*(M)B*(M)> = F B {k,t) (s» - (4) 

| Kim and Hughes (1997) considered values of magnetic Reynolds number Rm up to 10 5 . We use 
Rm up to 10 12 , restricted only by a limited numerical precision ~ 1CP 16 of floating-point numbers of 
double type in our computer code. 
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where F B (k,t) is the magnetic energy spectral function, (|B(k, t)| 2 ) = 2F B (k,t), and 
Hs(k, t) is the spectral function of the electric current helicity, (_B l *(k, t)ie l i l ti B l (k, t)) = 
HB(k,t). Functions -Fb(/c, t) and HB{k,t) can be obtained from functions Mi{x,t) and 
K(x,t), and vice verse, by using the three-dimensional Fourier transforms (Monin and 
Yaglom 1971). The problem is then to find the correlation function ([3]) of the magnetic 
field, or, alternatively, its Fourier version (j4]). 

The growing magnetic field B(x, t), amplified by dynamo action, satisfies the 
induction equation 

d t B = V x (v x B) + r/V 2 B, (5) 

where r\ is the magnetic diffusivity. Suppose that the velocity field (JTJ) (J2j) is given, 
i.e. kinetic energy kl{x) and kinetic helicity g(x) are given. It turns out that in this case, 
to find the properties of the growing magnetic field, one needs to solve two coupled linear 
homogeneous partial differential equations for functions Ml(x, t) and K(x,t) related to 
magnetic energy and magnetic helicity. These equations were first derived by Vainshtein 
and Kichatinov (1986) in the framework of the Kazantsev-Kraichnan model with non- 
zero kinetic helicity. Recently, Boldyrev et al. (2005) established that the system of 
Vainshtein-Kichatinov equations possesses a self-adjoint structure, which is somewhat 
similar to a two-component quantum mechanical "spinor" form with imaginary time. We 
solve these two self-adjoint linear homogeneous partial differential equations numerically 
by the fourth-order Runge-Kutta integration method and by matching the numerical 
solution to the analytical asymptotic solutions at x — * and x — > oo, for details see 
Malyshkin and Boldyrev (2007). A general solution for the growing magnetic field is 
a linear superposition of exponentially growing eigenmode solutions. Each eigenmode 
solution grows in time as exp(At), where A is the eigenmode growth rate. The self- 
adjoint structure of Boldyrev et al. (2005) equations guarantees that all growth rates 
are real. We find that, in analogy with quantum mechanics, there are two types of 
field eigenmodes: bound (spatially localized) and unbound (spatially non- localized). 
First, for growth rates A > A = g 2 (0)/[Kk(0) + 2rf\ the eigenf unctions are bound and 
correspond to "particles" trapped by the potential provided by velocity fluctuations. 
The bound eigenmodes have discrete growth rates, i.e., A = A n > Ao, n = 1,2,3.... 
As x — > oo the bound eigenfunctions decay exponentially to zero. Second, for A < A 
the eigenfunctions are unbound and correspond to "traveling particles". The unbound 
eigenmodes have a continuous spectrum of their growth rates, < A < Ao- The 
unbound eigenfunctions asymptotically become a mixture of cosine and sine standing 
waves as x — > oo. Eigenvalue Ao corresponds to the fastest growing unbound eigenmode. 
The properties of the magnetic field amplified by kinematic dynamo action are fully 
determined by all growing eigenmodes of the field. In particular, the magnetic energy 
spectral function Fs{k,t) is a linear sum of all energy spectral eigenfunctions, and the 
same is true for the electric current helicity spectral function -£/#(&, t). Thus, to find the 
properties of the growing magnetic field, it is sufficient to find growth rates and spectral 
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Figure 1. Growth rates X n of the bound magnetic eigenmodes, and Ao of the fastest 
growing unbound eigenmode. Plots A and B are for h = 1 and 0.01 respectively, while 
r] = 10 -6 (Rm ~ 10 6 ). Plots C and D are for h = 1 and 0.01 respectively, while 
j] = 10~ 12 (Rm - 10 12 ). 



eigenfunctions of all field eigenmodes 
3. Results and conclusions 

In this work we choose Gaussian velocity correlation tensor (J2J), 



k l = e , g 



where h is the kinetic helicity parameter, for which the condition — 1 < h < 1 must be 
satisfied jjj] Without loss of generality, the choice (J6j) means that all velocity fluctuations 
are of order unity, v ~ 1, and are on the scale of order unity, l ~ 1. Thus, the Reynolds 
number is of order unity and the turbulence has single-scale turbulent eddies. 

We study two cases for the magnetic Reynolds number, Rm. In the first case, 
Rm ~ 10 6 , which is achieved by choosing magnetic diflusivity rj = 10~ 6 . In the second 
case, Rm ~ 10 12 , which corresponds to our choice rj = 10~ 12 . These two cases for the 
Rm value are considered in combination with two cases for the kinetic helicity: first, 
a case when h — 1 in equation ([6]) and the kinetic helicity is large, and, second, a 
case when h = 0.01 and the kinetic helicity is small. Thus, in total we consider four 
cases for our choices of magnetic Reynolds number Rm and kinetic helicity parameter 
h. The resulting growth rates X n of the bound (localized) eigenmodes and Ao of the 
fastest growing unbound (non-localized) eigenmode are shown in figure HJ The growth 
rates are measured in the units of a turbulent eddy turn-over rate ~ vq/Iq ~ 1. For the 
cases Rm ~ 10 12 and h — 1, 0.01, the logarithmic-scale plots of the absolute values of 

§ While the energy spectrum FB(k,t) = (|B(k, t)| 2 )/2 is positive, the energy spectral eigenfunctions 
may be negative. 

Functions kl(x) and g{x) can not be chosen arbitrarily, their Fourier images G{k) and F(k) must 
satisfy the realizability condition |G(fc)| < F{k)/k (Moffatt 1978), which results in — 1 < h < 1. 
Analogously, functions Ml{x,€) and K(x,t) are restricted by the condition |iJs(fc,t)| < 2kFB(k,t). 
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Figure 2. Absolute values of magnetic energy spectral eigenfunctions for four selected 
bound eigenmodes (shown by the dotted, dash-dot, dashed and smooth solid lines), 
and for the fastest growing unbound eigenmode (shown by the red jagged spiky solid 
lines). All plots are for r\ = 10~ 12 (Rm ~ 10 12 ). The left and right plots are for the 
cases h = 1 and h — 0.01 respectively. 




magnetic energy spectral eigenfunctions are given in figure [2] The absolute values of 
the corresponding current helicity spectral eigenfunctions are plotted in figure [31 

In the case Rm ~ 10 6 and h = 1 there exist four bound magnetic field eigenmodes. 
They have growth rates Ai ~ 0.9991, A2 — 1.669, A3 ~ 2.681 and A4 ~ 3.443, which are 
shown on the plot A in figure [TJ The fastest growing unbound eigenmode grows at a 
rate Ao — 0.9943. In the case Rm ~ 10 6 and h = 0.01 the number of bound eigenmodes 
is again four, Ai ~ 0.1750, A 2 ~ 1.507, A 3 ~ 2.629 and A 4 ~ 3.430, while the fastest 
growing unbound eigenmode is Ao — 9.943 x 10~ 5 , refer to the plot B in figure [Q Next, 
in the case Rm ~ 10 12 and h = 1 there are seven bound eigenmodes, all shown on the 
plot C in figure [H Among these we select four bound modes Ai ~ 1.192, A 3 ~ 2.322, 
A6 — 3.499, A7 ~ 3.686, and we plot their spectral eigenfunctions by the dotted, dash- 
dot, dashed and smooth solid lines respectively on the left plots in figures [2] and [31 The 
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spectral eigenfunctions of the fastest unbound mode A — 0.9943 are shown by the red 
jagged spiky lines on these left plots. Finally, in the case Rm ~ 10 12 and h = 0.01 there 
are eight bound eigenmodes, refer to plot D in figure [U The spectral eigenfunctions 
of four selected bound modes Ai ~ 0.2948, A 3 ~ 1.694, A 7 ~ 3.494 and A 8 ^ 3.685 
are shown by the dotted, dash-dot, dashed and smooth solid lines on the right plots in 
figures [2] and [3j The spectra of the fastest unbound mode A — 9.943 x 10~ 5 are again 
shown by the red jagged spiky lines. As an interesting result, note that in figures [2] and [3] 
the envelopes of all bound eigenfunctions - not only the envelope of the fastest growing 
eigenfunction - agree with the Kazantsev spectral slope k^ 1 ^ 2 (compare to Kulsrud and 
Anderson 1992). 

Using the results presented in figures [THSl we make the following conclusions. First, 
the value of kinetic helicity parameter h does not have much effect on the fastest growing 
bound eigenmodes (e.g., compare modes A5, A6, A7 on plot C to modes A6, A7, As on 
plot D in figured]), but it does influence the unbound eigenmodes. This result is expected 
because helicity of velocity fluctuations does not propagate to small scales (Kulsrud and 
Anderson 1992), however, A = g 2 (0)/[fi; £ (0) + 2rf\ ~ h 2 v /l . Second, at a fixed kinetic 
helicity value, when the magnetic Reynolds number increases, the number of bound 
modes increases (refer to figure [[]). This is also expected because, the diffusion of 
magnetic field decreases with an increase of the magnetic Reynolds number. Third, let 
us refer to the spectra of the fastest growing unbound eigenmodes Ao, which are shown 
by the the red jagged spiky solid lines in figures [2] and [31 We see that when the value 
of the kinetic helicity drops by a factor of a hundred (from h = 1 to h = 0.01), the 
location of the peaks of these spectra shifts to larger scales by the same factor. Thus, 
the characteristic scale of eigenmode Ao is approximately equal to ~ lo/h, and this mode 
becomes large-scale when kinetic helicity is small, in agreement with the physical picture 
of large-scale dynamo action (Steenbeck et al. 1966). 

We thank Fausto Cattaneo for useful discussions. This work was supported by the 
NSF Center for Magnetic Self- Organization in Laboratory and Astrophysical Plasmas. 
SB is supported by the U.S. Dept. of Energy under the grant No. DE-FG02-07ER54932. 
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